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We examine cosmic microwave background constraints on variations in the recombination history 
of the universe. We use a very general extension to the standard model of recombination character- 
ized by two parameters, a and b, which multiply the overall rates of recombination and ionization 
and the binding energies of hydrogen respectively. We utilize WMAP temperature and TE cross- 
correlation data, along with temperature data at smaller scales from ACBAR and CBI, to place 
constraints on these parameters. The range of recombination histories which gives the best fit to 
the data is —0.6 < log(a) < 0.5 and 0.85 < b < 1.15, forming a diagonal region from low-a/high-6 
to high-a/low-b. We find Zrec = 1055 ± 25 and Az — 198 ± 5 at 68% confidence for the range of 
recombination models which are best able to reproduce the observed CMB power spectra. Standard 
recombination provides an acceptable fit to the data; while there is still room for non-standard 
recombination, the allowed variation is small. 



I. INTRODUCTION 

The cosmic microwave background (CMB) provides a 
snapshot of the early universe. The CMB radiation last 
interacted directly (i.e., non-gravitationally) with matter 
during the epoch of recombination, the period when the 
universe became cool enough for atomic nuclei to capture 
and hold electrons for the first time. Therefore, the CMB 
provides a sensitive test to any variations in the standard 
model of recombination. 

A number of authors have looked at CMB constraints 
on variations in standard recombination 0,0,0,^- Here 
we reexamine this question in the light of recent, high- 
quality data from the WMAP satellite ACBAR |||, 
and CBI JJ, following the parametrization presented in 
Ref. 0- The WMAP data represents a significant im- 
provement over data used in most previous calculations, 
and significantly tightens the resulting constraints. In ad- 
dition, we include the analysis of polarization data in the 
search for constraints on recombination. The only previ- 
ous analysis to use the recent WMAP data has been given 
by Bean, et al. who examined changes in the ioniza- 
tion history induced by Lyman-a photons and ionizing 
photons. Our approach is somewhat different; we have 
attempted a more general, model-independent investiga- 
tion, based on the discussion in Ref. jQj. Hence, the con- 
straints we derive are not tied to any particular model 
for altering recombination, but are potentially more gen- 
eral. Further, we have used additional data from ACBAR 
and CBI. In addition to using the TE cross-correlation 
data, we investigate the ability of the E-mode polariza- 
tion power spectrum to further constrain recombination. 

In Section ^] we describe a parameterization of the 
epoch of recombination and the details of how we ob- 
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tain constraints on the parameters. In Section IIIII we 
display the resulting constraints on recombination and 
discuss how the polarization power spectrum is affected 
by recombination. In Section llVI we present our conclu- 



II. METHODS 

A. Recombination Model 

The ionization fraction Xe is the fraction of ionized 
hydrogen present in the universe. 
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where Ue is the number density of free electrons and n 
is the number density of all hydrogen atoms, both ion- 
ized and neutral. The evolution of Xe is given by the 
ionization equation and describes the progress of recom- 
bination. The evolution of Xe is given by 
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Here C is Peebles' correction factor, R is the recombina- 
tion coefficient, (3 is the ionization coefficient, Bi and B2 
are the binding energies if the first two levels of hydro- 
gen, and Up is the number density of free protons plus 
hydrogen atoms. 

It is our goal to constrain recombination in a general 
way, rather than investigating specific mechanisms. To 
that end, we follow the framework presented in Ref. @ 
for modeling a general change in the evolution of the 
ionization fraction as a function of time. We modify the 
ionization history by inserting two parameters, a and 6, 
into Equation ||2Jl: 
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The constant a multiplies the overall rates of both re- 
combination and ionization, while the constant h multi- 
plies the binding energies of hydrogen. The values a = 1 
(or log(a) — 0) and 6=1 are the case of standard recom- 
bination. We neglect any corresponding change in helium 
recombination, as that effect is small and for most mod- 
els helium recombination is finished by the time hydrogen 
recombination begins. Very roughly, a change in a alone 
will change the duration of recombination, while keeping 
the onset of recombination fixed, while altering b shifts 
the onset of recombination to earlier or later redshifts. 



B. Data Sets and Priors 

Our goal is to find a range of recombination histories 
— that is, regions in a-b parameter space — which are ca- 
pable of predicting CMB power spectra compatible with 
the observed CMB power spectra. We create a grid in 
a-6, and for each point in this space, we allow a set of 
five cosmological parameters to vary, as well as the over- 
all normalization of the power spectrum. The vector of 
free parameters that we vary is z9={ilm, tJb, h, n^, trj, 
Qy. the matter density relative to critical 17™, the baryon 
matter density bjjj = f2f,/i^, the Hubble parameter h (in 
units of 100 km sec^^ Mpc^^), the spectral tilt rig, the 
optical depth to the surface of last scattering r^'/, and 
the overall normalization of the power spectrum Q. 

We assume priors for several of these cosmological pa- 
rameters. For the Hubble parameter, we impose the 
constraint h = 0.72 ± 0.08, consistent with estimates 
obtained from the Hubble Space Telescope (HST) Key 
Project and observations of Type la supernovae (SNe la) 
[Tol ITU . The baryon density is taken from BBN con- 
straints [III: = nbh"^ = 0.020 ± 0.002. Our estimate 
for the total matter density comes from supernovae stud- 
ies ^301 and x-ray observations of the gas mass fraction 
of relatively relaxed galaxies 0|. We combine these, us- 
ing generous error bars, to get VLm = 0.29 ± 0.10. These 
are all taken to be Gaussian, with uncertainties at the 
68% confidence level. For the spectral tilt we impose 
a uniform prior of 0.7 1.3. Values of beyond 

this range are disallowed. The reason for this is that Ug 
is generally not well known; indeed, the best estimates 
for the value of Ug come from the CMB itself [H [13 [H. 
Thus we use generous bounds which are theoretically mo- 
tivated. We do not include priors for the optical depth 
Tjii or the overall normalization of the power spectra Q. 
We additionally restrict our models to a fiat ACDM uni- 
verse (ritot = 1) and adiabatic initial density perturba- 
tions. 

We then use CMBFAST ligl| to produce a CMB 
anisotropy spectrum from our model, as a function of 
a, b, and 1}. For each point in a-b space we vary to 
minimize giving the combination of free parameters 
which is best able to reproduce the observed CMB power 
spectra for a given a and b. 

For our data we use a combined set consisting of the 



WMAP temperature and TE cross-correlation 
data, along with smaller scale temperature data from 
ACBAR 131 and CBI HIH. We follow the lead of Ref. 
23] in combining the WMAP temperature data with the 
higher £ data from ACBAR and CBI. A calibration uncer- 
tainty of 20% is assumed for ACBAR and 10% for CBI. 
Although ACBAR has a 3% beam width uncertainty, this 
effect is extremely small compared to the calibration un- 
certainty and we neglect it. We minimize using the 
combination of temperature and TE power spectra. 

III. RESULTS 

A. Effects of Varying the Recombination History 

Altering a or 6 from their standard values has a direct 
effect on how recombination proceeds. The easiest way 
to see this is by looking at the ionization fraction Xg as a 
function of time. The ionization fraction determines the 
optical depth, r due to Thomson scattering, 

T — — caTne{z){dt/dz)dz, (4) 
^0 

which in turn gives the visibility function, g{z): 

g{z) = e^'^dr/dz, (5) 

Figure ^ shows the ionization fraction versus redshift 
for a variety of recombination histories (a, b). In each 
case, the free parameters i9 are fixed at the values {ilm — 
0.29, uJb = 0.020, h = 0.72, = 1.0, r = 0.15}. As a 
increases, recombination has a longer duration, though 
it starts at about the same time. Beyond a certain limit, 
roughly log(a) = 1.5, recombination proceeds as rapidly 
as possible, and increasing a has no further effect because 
in the limit of high a, recombination proceeds in equilib- 
rium. As b increases, the primary effect is that the onset 
of recombination shifts to earlier epochs. A secondary ef- 
fect is that the duration of recombination increases with 
increasing b. 

The resulting visibility functions of the models used in 
Figure n are shown in Figure El As a decreases, the red- 
shift of recombination, Zrec, moves to more recent times, 
and the width Az of the last scattering surface (LSS) in- 
creases. (We define Zrec to be the redshift at which the 
visibility function is peaked, while Az is the full-width 
half-maximum of the visibility function) . The degeneracy 
of models with sufficiently high values of a is clearly seen. 
As b decreases, the peak of the visibility function and 
Zrec move toward lower redshifts, just as with decreasing 
a. However, with decreasing 6, unlike decreasing a, the 
width Az decreases. Thus, at the expense of allowing Az 
to change, changing a and b in opposite directions allows 
Zrec to remain relatively unchanged. 

Changes to the recombination history have a direct ef- 
fect on the temperature and polarization anisotropy spec- 
tra of the CMB. Decreasing a effectively increases the 
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FIG. 1: The ionization fraction, Xe{z), as a function of red- 
shift, z, for varying a, fixed b (top) and for varying b, fixed 
a (bottom), for Qm = 0.29, ujb = 0.020, h = 0.72, n« = 1.0, 
TRJ = 0.15. 

width of the last scattering surface Az, and to a lesser 
extent also decreases Zrec- The increase in Az leads to 
increased diffusion damping of the anisotropics at large 
scales so that the peaks in the temperature, polarization, 
and TE cross-correlation power spectra at large scales 
are all less pronounced. The decrease in z^ec shifts the 
features of all three power spectra to larger scales. 

Decreasing b primarily decreases z^ec and also leads to 
a slight decrease in the value of Az. The decrease in Zrec 
will again shift the features of the power spectra to larger 
scales. In this case the relative amplitudes of the peaks 
vary considerably, whereas when a is varied, the relative 
amplitudes remain fairly constant: as b is increased, the 
height of the even peaks relative to the odd peaks is also 
increased. This is due to the change in the ratio of pb to 
p-y caused by later recombination. Because Az increases 
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FIG. 2: As Fig. 1, for the visibility function g{z). 

only slightly as b is increased, diffusion damping at small 
scales for higher values of b is less significant. 



B. Allowed Region in a-b Space 

The region in a-b parameter space which is consistant 
with the CMB observations is displayed in Figure|2| This 
is the main result of our paper. The single best fit value 
is log(a) — —0.2, b = 1.0, and standard recombination 
lies well inside the la region. The allowed region of a-b 
space is considerably reduced compared to constraints 
using pre-WMAP data 01 ■ The 3a region extends only 
from -0.8 < log(a) < 1.2 and 0.70 <b< 1.25. The la 
regions are —0.6 < log(a) < 0.5 and 0.85 < 6 < 1.15. 
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FIG. 3: Allowed region of a-b space. Solid, long-dashed, and 
short-dashed curves enclose the la, 2a, and 3a regions. 



C. Characteristics of the Allowed Regions 

It is interesting to examine the manner in which the 
various cosmological parameters change in order to com- 
pensate for changes in a and b. 

The shape of the CMB power spectrum is sensitive to 
changes in each of the cosmological parameters In- 
creasing the baryon density suppresses the even peaks 
with respect to the odd peaks because the increased 
baryon mass contributes more gravity to the potential 
wells, lessening the effect of the rarefaction phase of the 
acoustic oscillations. Increasing the Hubble constant will 
lead to matter-radiation equality happening earlier and 
an increased expansion rate, thus the power spectrum 
will tend to have suppressed even peaks and move to 
slightly lower values of £. An increase in the total matter 
density keeping the total density fixed at SI tot — 1, or 
equivalently a decrease in the dark energy density, again 
moves matter-radiation equality to earlier times, which 
tends to move the spectra to lower £ and increase the late 
ISW effect. An increase in the scalar spectral index gives 
more power to smaller scales, and will tilt the CMB power 
spectrum, increasing the amplitude at large £ compared 
to small £. Finally, increasing the optical depth due to 
rcionization will decrease the power at small scales by a 
factor e^^'^'^' . 

Figure 0] shows how the best-fit cosmological parame- 
ters change for all points falling inside the Su boundary. 
In producing this figure, we first determine the mean and 
1(7 variation for these parameters within the la region 
in a — 6 space. Then we calculate the best-fit value of 
each parameter for a given value of a and b. The open 
squares with x 's show points for which the best-fit value 
for a given cosmological parameter lies within ilcr of 
the mean value. The solid squares show points for which 
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FIG. 4: How the best-fit cosmological parameters change with 
a and b for recombination histories falling inside the 3a re- 
gion. The squares with an x have "average" values for a 
given parameter, the solid squares have "high" values, and 
the open squares have "low" values. (See text for the precise 
definition) . 



the best-fit value of the given parameter is higher than 
this (i.e., more than la above the mean), and the open 
squares show points for which the best-fit value is lower 
than la below the mean. 

The baryon density uJb shows a strong tendency to in- 
crease both as a is increased and as b is increased. Simi- 
larly, the Hubble constant h increases both with increas- 
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ing a and increasing b. The scalar spectral index in- 
creases with decreasing a but there does not appear to 
be any correlation of Ug with changing b. The optical 
depth Tjii appears to increase with decreasing 6, but this 
effect is not nearly so monotonic as the correlations in 
the previous three parameters. There does not appear 
to be any significant correlation of tri with changes in 
a. For firm there does not appear to be much correlation 
with changes in a or b. 

To summarize these effects, an increase in a is com- 
pensated by corresponding increases in uji, and h and a 
corresponding decrease in Us- An increase in b is com- 
pensated by increases in LUh and h and possibly a decrease 
in TRi. 

The corresponding changes in the cosmological param- 
eters as a and b are changed are due to competing effects. 
As a is increased, the peaks at high £ are less damped, be- 
cause the width of recombination is decreased, and the 
peaks are shifted to higher £ because recombination is 
moved to higher redshift. The decrease in Ug lowers the 
amplitude of the spectrum at small scales, counteracting 
the increased power at small scales due to a narrower sur- 
face of last scattering. An increase in h tends to move the 
peaks to larger scales, counteracting the effect of a higher 
Zrec- Similarly, increasing b will enhance the height of the 
even peaks by increasing Zrec, thus lowering the ratio of 
the baryon to photon density. This is compensated by 
a higher values of h and wt, which tend to suppress the 
height of the even peaks. The increased h will again move 
the peaks to slightly larger scales, even as an increase in 
b also increases Zrec and shifts the spectrum to smaller 
scales. 

It is also instructive to look at the evolution of the ion- 
ization fraction function of the redshift z for those 
recombination histories which provide a good fit to the 
data. Figure |5l shows Xe versus redshift for a selection of 
recombination histories in the la, 2a, and 3a regions. In 
each case, the cosmological parameters used in calculat- 
ing Xe{z) were those which produced the best fit to the 
CMB data. 

From Figure \E\ it can be seen that the common fea- 
ture in all recombination histories capable of producing 
the observed CMB is their ending point. The recombina- 
tion histories start at different times and have different 
durations, but all end at approximately the same era. 
For the la region, the average redshift at which the ion- 
ization fraction falls to 1% is z{xe — 0.01) = 856 ± 20. 
By contrast, the point where the ionization fraction falls 
to 50% has a much larger standard deviation: z(xe — 
0.50) = 1237 ± 44. All error estimates are at the 68% 
confidence level. 

We find a similar but smaller effect looking at the visi- 
bility function. Figure [S] shows the corresponding visibil- 
ity functions g{z) for the la, 2a and 3a regions. For the 
la region, we find that the average redshift of the peak 
of the visibility function is Zrec = 1055 ± 25, while the 
average width is Az = 198 ± 5. If we look at the point 
where the visibility function falls to 10% of its peak value. 
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FIG. 5: Ionization fraction Xe{z) vs. redshift z (top) and visi- 
bility function g{z) vs. redshift z (bottom) for recombination 
histories falling within the la, 2a, and 3a regions of the anal- 
ysis. The thicker dashed line is standard recombination, and 
it falls within the la region. It is shown on the other two 
graphs for comparison purposes. 

we again find a somewhat tighter constraint than at the 
peak, z{g/g,nax = 0.10) = 811±I9, though the difference 
is not quite as great as that found using Xe- 

The shapes of the visibility functions stray further and 
further from that of standard recombination in the 2a 
and 3a regions. Although some curves in the 3a re- 
gion have a similar Zrec to standard recombination, these 
curves tend to have a significantly different Az. This 
suggests that we are not placing a constraint on z^ec or 
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Az individually, but jointly. 
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FIG. 6: Redshift of recombination Zrec vs. width of the 
last scattering surface for recombination histories falling 
within the Icr, 2cr, and Sa regions of the WMAP+ TT+TE 
analysis. The apparent lines in this graph are due to the finite 
grid spacing in a and b. 

In Figure IHI we have taken points from the Icr, 2cr, 
and 3(7 regions of the allowed a-h region, calculated their 
visibility functions, and plotted the resulting values of 
Zrec versus Az. Although we are not drawing confidence 
regions in this z^ec — Az space, it can be clearly seen 
that the recombination histories which are best able to 
reproduce observations lie in a cluster, with successively 
poorer histories lying further out. As noted earlier, for 
all of the ionization histories inside the la region, z^ec = 
1055 ± 25, and Az = 198 ± 5. 

In Figure we show the temperature, TE cross- 
correlation, and E-mode polarization power spectra for 
recombination histories falling within the Icr region. The 
positions of the peaks in all three plots seem very well 
constrained, but the heights of the peaks differ somewhat. 
The Cj"^ and Cj^ power spectra match the data very 
well, which is unsurprising since that was our criteria. 
Given the likely precision of data on the Cf^ power spec- 
trum in the near future, we conclude that it is unlikely 
that measurements of the E-mode polarization would be 
able to break any degeneracies in the allowed a-h param- 
eter space and significantly improve the constraints on 
recombination. 
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FIG. 7: Cf^, Cf ^ and Cf ^ vs. I for recombination histories 
falling inside the la region. The thicker dashed line in each 
plot is standard recombination. 



IV. CONCLUSIONS 

We have examined the new data from the WMAP mis- 
sion, along with data from higher i experiments, ACBAR 
and CBI, in an attempt to constrain recombination. We 
have included analyses of both the temperature power 
spectrum Cj^ and the TE cross-correlation power spec- 
trum Cj^ . In the general parameterization that we used, 
the Icr allowed region in a-h space was significantly con- 
tracted compared to earlier attempts with less precise 
data (jSi]). 



We find that the allowed region falls within —0.6 < 
log(a) < 0.5 and 0.85 < 5 < 1.15. The average width and 
position of the surface of last scattering for recombination 
histories falling within the Icr region are z^ec = 1055 ± 25 
and Az = 198 ±5. Note that the WMAP team estimated 
values for Zrec and Az, namely [H zJ^^ap ^ loss^^ and 
Az^MAP ^ 294 ± 2. Our values for Zrec and Az do not 
mean the same thing as those reported by the WMAP 
team. In the latter case, standard recombination was as- 
sumed, and these values simply give the best-fit recom- 
bination history obtained by allowing the cosmological 
parameters to vary. In our case, we allow variations in 
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the recombination history itself, so we expect a larger 
range of values for z^ec a-nd Az. In fact, our values for 
Zrec and Az are quite similar to the WMAP values; this 
is a reflection of the fact that standard recombination is 
an excellent fit to the data and is approximately in the 
center of the distribution of the allowed variant recombi- 
nation histories. 

The changes in a and h are largely compensated by 
changes in other parameters. The largest effects are that 
an increase in a is matched by corresponding increases 
in Wft and h and a corresponding decrease in n^, while 
an increase in b is compensated by increases in ujh and 
h. Tighter constraints on these parameters, especially h, 
would allow us to narrow the range of allowed recombi- 
nation histories. 

The result of these competing effects — changes in a 
and h being partially compensated for by changes in the 
free cosmological parameters — is that the anisotropy 
power spectra from the la region all share very simi- 
lar features. In general, the positions of the peaks are 
strongly constrained, while the amplitudes of the peaks 
are less so. The range in E-mode polarization power spec- 
tra corresponding to recombination histories from this 
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